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Transformación  de  una  función  de  entrada  fe  en 
otra  de  salida  fs  por  la  acción  de  un  sistema  X. 


Introducción 


El  análisis  de  las  imágenes  de  satélite  puede  abor¬ 
darse  desde  un  punto  de  vista  puramente  visual,  con 
los  mismos  criterios  que  se  aplican  a  la  fotointerpre- 
tación,  aunque  con  los  matices  que  introduce  el  ha¬ 
ber  sido  adquiridas  desde  plataformas  orbitales  en  lu¬ 
gar  desde  aeronaves  volando  a  baja  altura.  Sin  em¬ 
bargo,  el  potencial  de  información  que  proporcionan 
las  imágenes  de  satélite  requiere  un  estudio  más  pro¬ 
fundo  que  solo  puede  realizarse  mediante  análisis  di¬ 
gital.  Para  ello,  las  imágenes  deben  ser  consideradas 
objetos  matemáticos  susceptibles  de  ser  manipuladas 


Sistemas 


Un  sistema,  es  una  estructura,  dispositivo,  esquema 
o  procedimiento  -real  o  abstracto-  que  relaciona,  en 
un  tiempo  de  referencia  dado,  una  entrada,  causa  o 
estímulo  de  materia,  energía  o  información  y  una  sa¬ 
lida,  efecto  o  respuesta  de  información,  energía  o  ma¬ 
teria. 

Esta  definición  genérica,  destinada  en  su  origen  al 
campo  de  la  hidrología,  engloba  la  práctica  totalidad 
de  los  espacios  de  aplicación  de  la  ingeniería  de  sis¬ 
temas  y,  en  particular,  al  análisis  digital  de  imágenes. 
En  el  ámbito  de  la  teledetección,  el  factor  tiempo  in¬ 
cluido  en  la  definición  suele  carecer  de  importancia 
al  aplicarse  generalmente  a  imágenes  estáticas.  Sin 
embargo,  en  otros  campos,  sí  que  es  algo  para  tener 
en  cuenta,  como  sucede  en  el  análisis  y  el  manejo  de 
imágenes  en  movimiento,  como  por  ejemplo  en  visión 
artificial  o  en  televisión. 


en  sistemas  informáticos  para  poder  realizar  sobre 
ellas  operaciones  matemáticas,  a  partir  de  las  cuales 
se  puedan  deducir  propiedades  inherentes  a  la  ima¬ 
gen  que  en  el  simple  análisis  visual  pasan  inadverti¬ 
das. 

Para  abordar  el  análisis  digital  de  las  imágenes  es 
muy  útil  considerar  el  conjunto  [imagen  bruta  -  sis¬ 
tema  de  tratamiento  -  imagen  derivada]  como  un  sis¬ 
tema.  En  el  apartado  siguiente  se  realizará  una  apro¬ 
ximación  a  la  ingeniería  de  sistemas. 


Es  importante  considerar  que  el  sistema  puede  estar 
constituido  por  uno  o  más  elementos  reales  susten¬ 
tado  en  soportes  físicos,  o  bien  ser  una  estructura 
abstracta  que  transforma  entradas  en  salidas,  y  que 
éstas  pueden  ser  de  la  misma  naturaleza  o  no.  En  el 
ámbito  de  las  finanzas,  una  entrada  de  naturaleza  in¬ 
formativa,  que  puede  considerarse  abstracta,  puede 
producir,  tras  las  operaciones  bursátiles  necesarias, 
una  salida  basada  en  el  pánico  que  resulte  en  una 
caída  masiva  de  las  cotizaciones,  aspecto  que  sí  tiene 
naturaleza  real  y  contable.  En  otras  ocasiones,  la  en¬ 
trada  y  la  salida  son  del  mismo  tipo,  como  sucede  en 
el  campo  de  los  estudios  hidrológicos.  La  entrada  a 
una  cuenca  hidrográfica  puede  ser  la  precipitación 
registrada  durante  un  aguacero  o  en  un  episodio  de 
lluvia  extrema,  y  la  salida  es  el  caudal  aforado  o  la 
cota  de  la  lámina  de  agua  en  la  llanura  de  inundación 
medida  en  un  punto  sensible  de  la  cuenca  aguas 
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abajo.  Bajo  la  teoría  de  sistemas  se  analizan  aspectos 
muy  diversos  de  la  vida  cotidiana,  desde  la  sociología 
pasando  por  la  geografía,  la  economía  y  la  informa¬ 
ción  hasta  numerosas  especialidades  de  la  ingeniería. 

En  términos  generales,  un  sistema  está  constituido 
por  tres  componentes: 

-i  La  entrada  X(t). 

— i  El  operador  o  función  de  transferencia  X 
-i  La  salida  Y(t). 

La  relación  existente  entre  la  entrada  y  la  salida  se 
expresa  mediante: 

Y(t)=x{x(t)}  (1) 

que  indica  que  Y{t)  es  la  respuesta  de  X  ante  la  en¬ 
trada  X{t). 

Los  sistemas  pueden  ser  lineales  o  no  lineales.  Aun¬ 
que  en  el  ámbito  del  tratamiento  digital  de  imágenes 
existen  aplicaciones  basadas  en  estos  últimos,  en  la 
mayoría  de  los  casos,  los  que  se  utilizan  son  los  linea¬ 
les.  Para  que  un  sistema  sea  lineal  debe  cumplir  dos 
propiedades: 

-i  Principio  de  proporcionalidad:  si  a  la  entrada  Xi  al 
sistema  le  corresponde  la  salida  Yi,  la  entrada 
C-Xi  debe  producir  una  salida  C-Yi,  siendo  C  una 
constante. 

-i  Principio  de  aditividad:  si  dos  entradas  individua¬ 
les  Xi  y  X2  producen,  respectivamente,  dos  salidas 
Yi  e  Y 2,  la  entrada  X1+X2  debe  producir  la  salida 
Y1+Y2. 

Estas  dos  propiedades  se  pueden  combinar  escri¬ 
biendo: 


Carlos  Pinilla 


X  (Ci  Xj  +C2  x2)  = 

=  X  (Ci  Xi)  +  X  (C2  X2)  =  C,  X  (Xi)  +  C2  X  (X2)  = 

=  C1Y1 +C2Y2  (2) 

que  significa  que  el  resultado  de  aplicar  el  sistema  a 
la  suma  de  dos  entradas  amplificadas  es  el  mismo  que 
si  se  suman  los  resultados  amplificados  de  las  salidas 
de  cada  una  de  las  entradas  sometidas  al  sistema.  En 
palabras  más  sencillas,  considerando  el  filtrado  de 
imágenes  como  un  sistema  lineal  y  teniendo  en 
cuenta  solamente  la  propiedad  de  aditividad,  se 
puede  decir  que  el  filtrado  de  la  suma  de  dos  imáge¬ 
nes  produce  el  mismo  resultado  que  la  suma  de  los 
filtrados  de  cada  una  de  ellas.  Más  adelante  se  verá 
que  también  hay  filtros,  así  como  otras  operaciones 
interesantes  en  teledetección,  que  no  son  lineales. 

Los  instrumentos  de  adquisición  de  imágenes  son, 
por  lo  general,  sistemas  lineales  invariables  en  el 
tiempo,  en  los  cuales  la  magnitud  de  la  salida  no  de¬ 
pende  del  tiempo  sino  solo  de  la  entrada.  En  con¬ 
creto,  en  un  sistema  formador  de  imágenes,  la  en¬ 
trada  es  la  radiancia  de  la  escena,  el  sistema  está 
constituido  por  el  sensor  y  la  salida  es  la  irradiancia 
sobre  el  plano  focal  del  sensor  o,  si  además  de  la  óp¬ 
tica  se  considera  la  electrónica  del  instrumento,  su 
traducción  a  voltaje  o  niveles  digitales  de  la  imagen 
formada.  En  teledetección  se  considera  a  menudo  la 
atmósfera  como  una  parte  más  del  sistema,  ya  que  su 
presencia  modula  la  señal  procedente  de  la  superfi¬ 
cie. 
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La  imagen  digital 

En  su  definición  académica,  la  imagen  es  la  reproduc¬ 
ción  de  la  figura  de  un  objeto  por  la  combinación  de 
los  rayos  de  luz  que  proceden  de  él.  El  diccionario 
realmente  no  aclara  gran  cosa,  puesto  que  propor¬ 
ciona  una  descripción  circular  del  concepto  cuando 
habla  de  reproducir  la  figura.  Más  allá  de  las  defi¬ 
niciones  formales,  aquí  lo  que  interesa  es  su  des¬ 
cripción  física,  y  desde  el  punto  de  vista  óptico,  una 
imagen  es  una  distribución  bidimensional  de  energía 
electromagnética,  particularmente  circunscrita  al  es¬ 
pectro  visible.  Aun  así,  la  teledetección  amplía  la  no¬ 
ción  a  cualquier  otro  intervalo  de  longitudes  de  onda 
del  espectro  sin  ningún  problema,  y  de  este  modo, 
aunque  no  sea  visible  por  el  ojo  humano,  se  puede 
configurar  perfectamente  una  imagen  a  partir  de  la 
radiancia  espectral  extendida  a  cualquier  tipo  de  ra¬ 
diación  electromagnética.  De  manera  más  sucinta,  la 
imagen  digital  es  un  modelo  digital  de  radiancia. 

Considérese  una  fuente  radiante  cuya  distribución  es¬ 
pectral  de  radiación  esté  dada  por  una  función  o,.  Si 
la  energía  que  procede  de  ella  incide  sobre  una  su¬ 
perficie  de  reflectancia  espectral  p*,,  la  energía  de¬ 
vuelta  para  cada  longitud  de  onda  será: 

C.  =  cxPa.  (3) 

y  la  devuelta  para  el  intervalo  de  longitudes  de  onda 
comprendidos  entre  Xi  y  X2: 

eM2]=£\p  (4) 


Esto  quiere  decir  que  la  energía  que  puede  regis¬ 
trarse  como  consecuencia  de  la  reflexión  de  la  radia¬ 
ción  sobre  una  superficie  bajo  observación  depende 
de  la  composición  de  la  luz  incidente,  de  la  reflectan¬ 
cia  espectral  de  la  superficie  y  del  intervalo  de  longi¬ 
tudes  de  onda  estudiado.  Pero  es  importante  consi¬ 
derar,  además,  que  esta  energía  -y  en  definitiva,  la 
imagen  formada  a  partir  de  su  registro-  dependen  de 
cada  posición  en  el  terreno,  puesto  que  cada  punto 
puede  tener  una  reflectancia  diferente.  De  este  modo 
se  puede  proponer  un  modelo  matemático  para  una 
imagen  continua,  monocroma  y  fija,  esto  es,  conside¬ 
rando  una  sola  longitud  de  onda  y  además  prescin¬ 
diendo  de  cambios  o  de  movimiento: 

z  =  f{x,y )  (5) 

siendo  x,y,z  E  9Í,  esto  es,  pudiendo  tomar  cual¬ 
quiera  de  las  variables  valores  reales  arbitrariamente 
próximos. 

Es  posible  normalizar  la  imagen  adoptando  un  rango 
de  existencia  para  la  variable  z  de  modo  que  0  <  z  < 
1  y  atribuyendo  el  color  negro  al  0  y  el  blanco  al  1. 
Entre  ambos  extremos  se  distribuirá  toda  la  gama  de 
grises  proporcional  al  valor  de  la  imagen  en  cada 
punto. 

Un  sistema  de  tratamiento  digital  de  imágenes  no 
puede  procesar  una  imagen  de  este  tipo,  puesto  que 
las  posiciones  de  sus  puntos  difieren  infinitesimal¬ 
mente  en  distancia  y  su  tamaño  en  consecuencia 
tiene  al  infinito.  Se  hace  necesario  discretizar  la  ima¬ 
gen.  Se  entiende  por  imagen  discreta  aquella  cuyas 


Teledetección:  teoría  y  aplicaciones 


06-4 


La  imagen  digital 


posiciones  planimétricas  son  discretas,  es  decir,  que 
solo  pueden  ocupar  ciertas  posiciones  y  no  las  inter¬ 
medias.  Para  construir  una  imagen  así  a  partir  de  una 
real  es  preciso  muestrear  el  espacio  objeto  generando 
una  malla  regular  y  atribuir  un  solo  valor  de  brillo  a 
cada  nodo,  promediando  la  iluminación  procedente 
de  la  región  dominada  por  cada  uno  de  esos  puntos. 
El  modelo  matemático  para  esta  imagen  será: 

z  =  f(c,l )  (6) 

siendo  ahora  c,l  E  Z  y  z  E  %  o  sea,  las  posiciones 
sobre  el  plano  son  números  enteros,  que  pueden  or¬ 
ganizarse  en  columnas  (c)  y  líneas  (/),  aunque  el  brillo 
de  cada  punto  todavía  puede  adoptar  cualquier  valor 
real  entre  0  y  1.  Así  pues,  en  la  imagen  discreta,  las 
coordenadas  de  los  puntos  son  discretas  y  la  lumi- 
nancia  puede  variar  de  forma  continua.  Si  las  coorde¬ 
nadas  ( c , /)  cumplen  que  0  <  c  <  C- 1  y  0 < / < L- 1  en¬ 
tonces  la  imagen,  además  de  discreta,  es  finita  y 
puede  representarse  por  una  matriz  de  dimensiones 
CxL,  siendo  C  el  número  total  de  columnas  y  L  el  de 
líneas  de  la  imagen.  Cada  celda  de  la  matriz  corres¬ 
ponderá  a  una  celda  de  imagen.  Pero  sigue  siendo 
esta  imagen  imposible  de  someter  a  tratamiento  digi¬ 
tal  debido  a  la  continuidad  de  los  valores  de  luminan- 
cia.  Por  ello  es  preciso  hacer  discretos  también  los 
valores  de  luminancia.  De  esta  manera  se  enuncia  el 
modelo  de  imagen  digital: 

ND  =  f(c,l )  (7) 

donde  ND  significa  Nivel  Digital,  que  es  el  brillo  o 
valor  de  gris  discreto  adoptado  en  sustitución  de  la 
variable  continua  anterior  z.  Ahora,  en  consecuencia, 
c,  l,  ND  e  Z .  Finalmente,  para  que  un  sistema  infor¬ 
mático  pueda  tratar  estas  imágenes  es  necesario  im¬ 
poner  una  condición  más:  el  nivel  digital  también 


debe  adoptar  valores  enteros  en  la  forma: 

0  <  ND  <  2k  -  1  (8) 

y  por  lo  tanto  tener  un  valor  mínimo  igual  a  0  y  unos 
valores  máximos  que  dependerán  de  la  resolución  ra- 
diométrica  del  instrumento  de  adquisición,  tal  como 
muestra  la  Tabla  1. 


Tabla  1.  Número  de  ND  en  función  del  número  de  bits. 


Núm. 

Núm. 

Núm. 

Núm. 

bits 

ND 

bits 

ND 

1 

2 

9 

512 

2 

4 

10 

1024 

3 

8 

11 

2048 

4 

16 

12 

4096 

5 

32 

13 

8192 

6 

64 

14 

16384 

7 

128 

15 

32  768 

8 

256 

16 

65536 

Se  ha  dicho  anteriormente  que  la  imagen  de  satélite 
puede  explorar  regiones  de  longitudes  de  onda  del 
espectro  más  allá  del  visible.  Esto  hace  que  frecuen¬ 
temente  las  imágenes  almacenen  la  información  ob¬ 
tenida  del  terreno  en  distintas  bandas,  cada  una  de 
ellas  reservada  a  un  intervalo  espectral  concreto.  En 
consecuencia,  la  imagen  en  realidad  no  contiene  una 
sola  capa  de  información,  sino  tantas  como  bandas. 
Se  habla  entonces  de  una  imagen  multiespectral. 
Cada  celda  conserva,  pues,  tantos  niveles  digitales 
como  bandas  tenga  la  imagen.  El  modelo  matemático 
para  la  imagen  multiespectral  será: 

ÑD  =  f(c,l)  (9) 

es  decir: 


Proceso  de  discretización  de  la  imagen 
e  imagen  discreta  resultante. 
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238 

05 

255 

221 

0 

68 

136 

17 

170 

119 

68 

221 

0 

238 

136 

0 

255 

119 

255 

65 

170 

136 

238 

238 

17 

221 

68 

119 

255 

86 

170 

119 

221 

17 

136 

La  imagen  digital. 


La  imagen  multiespectral. 


NDl  ' 

nd2 

= 

El  conjunto  de  niveles  digitales  de  cada  celda  se  de¬ 
nomina  vector  de  características  y  representa  una 
posición  en  un  espacio  espectral  de  tantas  dimensio¬ 
nes  como  bandas  tenga  la  imagen  multiespectral. 


La  función  de  dispersión  puntual 


En  un  sistema  formador  de  imágenes,  el  nivel  digital 
atribuido  a  una  celda  concreta  corresponde  al  prome¬ 
dio  de  los  valores  de  radiancia  de  la  región  del  te¬ 
rreno  homologa.  Se  puede  hablar  de  espacio  objeto 
y  espacio  imagen  para  distinguir  entre  el  terreno  y  la 
imagen  formada  en  el  plano  focal  del  instrumento  de 
adquisición.  Desde  la  óptica  de  la  ingeniería  de  siste¬ 
mas  puede  considerarse  el  sistema  de  adquisición  de 
imágenes  como  el  conjunto  formado  por  el  espacio 
objeto,  la  atmósfera  interpuesta  entre  el  terreno  y  el 
sensor  y  el  propio  instrumento  formador  de  imáge¬ 
nes. 

Abandonando  por  un  momento  la  formación  de  imá¬ 
genes  bidimensionales  y  recurriendo  a  las  formula¬ 
ciones  clásicas  de  teoría  de  señales,  se  define  un  im¬ 
pulso  unitario  o  función  delta  de  Dirac1  como  una 
función  pico  de  intensidad  o  amplitud  infinita  y  una 
duración  infinitesimal.  Puede  definirse  el  impulso 
unitario  como  el  límite  de  un  pulso  rectangular  de 
duración  s  y  amplitud  1/s  cuando  s— >-0: 


S(í)  =  lím  — 

£^°  s  (11) 

La  definición  matemática  del  impulso  unitario  viene 
dada  por: 


«W 


0  V  t  ^  t0 

oo  t  =  0 


00  8 

|  8(í)  dt  =  |  S(t)  dt  =  1 

^  — 00  —8 


Y  para  el  instante  to: 


8(í-f0)  =  0  \/t^t0 

oo 

1 8(t-r0)  dt  =  1 


8  >  O 


(12) 


(13) 


El  impulso  unitario  tiene  propiedades  muy  interesan- 


Definición  de  la  función  5  de  Dirac  como  límite  del 
pulso  unitario. 


1  El  Nobel  de  física  inglés  Paul  Dirac  descubrió  el  electrón  y  predijo  la  existen¬ 
cia  del  positrón. 
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tes  en  el  tratamiento  de  señal  en  general,  y  de  imáge¬ 
nes  digitales  en  particular.  El  producto  del  impulso 
por  una  función  de  prueba  cualquiera  equivale  a  apli¬ 
carlo  solo  en  el  instante  en  el  que  el  impulso  no  es 
nulo: 

5(?  -  0/(0  Ü/(Ü  (14) 

Y  como  consecuencia  de  ello, 

00  00 

J  5(í-O/(0  dt=  {  6(?-0/(0  dt  = 

(15) 

00 

=  f(* o){5(r_0  dt 


Es  interesante  considerar  que  el  impulso  unitario 
también  se  puede  definir  como  la  derivada  de  la  fun¬ 
ción  escalón  unitario  de  Heaviside : 

8(0  = 

dt  (17) 

definiéndose  tal  función  como  sigue: 


f  u(t)f{t)dt  =  f  f(t)dt 

o  (i8) 

Es  decir,  cuando  se  aplica  la  función  de  Heaviside  a 
una  señal  de  prueba  sólo  extraen  los  valores  positivos 
de  esta  última.  Para  ello,  debe  ser: 


o  lo  que  es  lo  mismo,  si  se  introduce  el  intervalo 


J  5(í_O/(0  dt=  J  5(t)/(t  +  0  dx  = 


—00 


(16) 

00 

=  { 5(0/(í+O  dt=f(s) 


Es  decir,  la  integración  del  producto  de  un  impulso 
en  el  entorno  del  instante  to  por  una  función  o  señal 
da  como  resultado  el  valor  de  la  función  en  ese  ins¬ 
tante.  Esta  peculiaridad  se  conoce  como  propiedad 
de  exploración  o  muestreo  del  impulso  unitario. 


í  1  Vi  >  O 
|o  Vi  <  o 


(19) 


razón  por  la  que  se  denomina  escalón  unitario.  Para 
un  instante  genérico  to: 

u(t-t0)=< 


La  función  de  Heaviside  sólo  extrae  los  valores  de 
cualquier  función  de  prueba /(i),  a  partir  del  instante 
considerado  y  no  está  definida  para  i  =  to. 


1  Vi  >  i0 
O  Vi  <  in 


Mediante  ella,  el  impulso  unitario  permite  extraer,  re¬ 
tener  o  cribar  el  valor  de  una  señal  en  el  instante 
considerado.  Por  ello  sirve  para  muestrear  ( sample ) 
el  comportamiento  de  una  función  extendida  durante 
un  intervalo.  Es  la  forma  de  convertir  en  digitales  y 
discretas  señales  analógicas  continuas. 


En  el  contexto  del  tratamiento  digital  de  imágenes,  el 
concepto  de  impulso  unitario  ha  de  extenderse  a  dos 
dimensiones,  y  además  será  necesario  enunciarlo  de 
forma  discreta  y  no  continua.  El  impulso  unitario  bi- 
dimensional  puede  asimilarse  a  un  punto  o  celda  del 
máximo  brillo  sobre  un  fondo  completamente  oscuro. 


ü{t) 

A 


1 


O 


> 

t 


Escalón  unitario  de  Heaviside. 
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Función  de  dispersión  puntual 
en  la  formación  de  la  imagen. 


Cualquier  perfil  radiométrico  que  contenga  a  esa 
celda  brillante  tendrá  un  valor  0  hasta  llegar  a  ella, 
en  cuyo  momento  alcanzará  el  máximo  posible  para 
volver  a  0  de  nuevo.  La  duración  infinitesimal  en  las 
nuevas  circunstancias  es  un  segmento  espacial  del  ta¬ 
maño  de  una  celda  y  la  magnitud  infinita  se  convierte 
ahora  en  nivel  de  brillo  máximo  del  sistema. 

Si  la  entrada  al  sistema  formador  de  imágenes  fuese 
un  impulso  unitario,  en  ausencia  de  elementos  per¬ 
turbadores  entre  el  espacio  objeto  y  el  espacio  ima¬ 
gen,  la  imagen  formada  debería  ser  también  un  im¬ 
pulso  unitario,  es  decir,  la  imagen  de  una  celda  de 
terreno  brillante  habría  de  ser  un  píxel  brillante.  En 
la  práctica  esto  no  sucede  así  por  varias  razones  en¬ 
tre  las  que  las  más  importantes  son  la  presencia  de 
atmósfera,  con  efecto  dispersivo  debido  a  los  aeroso¬ 
les  y  moléculas  que  la  constituyen,  y  la  óptica  del  ins¬ 
trumento  de  adquisición,  que  puede  introducir  abe¬ 
rraciones  y  otros  efectos  dispersivos.  El  resultado  es 
que  la  imagen  formada  recoge  principalmente  el  pixel 
brillante,  pero  también  forma  un  halo  a  su  alrededor 
de  mayor  o  menor  entidad  dependiendo  del  alcance 
de  la  dispersión  operante.  La  forma  en  la  que  se  de¬ 
grada  la  señal  unitaria  original  se  conoce  como  Fun¬ 
ción  de  Dispersión  Puntual  (Point  Spread  Function, 
PSF). 

Es  necesario  en  este  punto  introducir  otra  operación 
matemática  para  comprender  cómo  se  forman  real¬ 
mente  las  imágenes.  Se  trata  de  la  integral  de  convo- 
lución.  Se  define  la  operación  de  convolución  de  dos 
funciones  f\(t)  y  fi(t)  como  la  nueva  función /(í): 


/,(0=o| 

/,(')= oj 


Vi  <0 


(22) 


sería: 


f(t)  =  ¡  fi(x)f2(t-x)dx 
0 


(23) 


La  operación  de  convolución  tiene  las  propiedades 
conmutativa  y  asociativa.  Así  pues: 

/i  M ' ®  fi  M  =  fi  (0  ®  /i  (0  (24) 


y 

WO®  /2  (01®  fAt)= /i(0®  [/2(0®  fM  (25) 

Cuando  una  de  las  dos  funciones  es  un  impulso  uni¬ 
tario,  se  tiene: 


/(f)08(f)=  j  f(x)  d(t  - x)  dx  (26) 

—oo 

y  por  su  propiedad  conmutativa: 

/(í)®8(f)  =  8(í)®/(f)  (27) 

Echando  mano  de  la  propiedad  de  muestreo  del  im¬ 
pulso  unitario  vista  en  (16): 


f{t)=f\{t)®f2{t)=\í\{x)f2{t-x)dx  (21)  j  5(7)/’(r  +  t0 )  dt  =  f(t0 ) 


que  en  el  caso  particular  de  ser 


se  tiene  que 


(28) 
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8  (t)®f(t)=  j8(x)f(t-x)dx  =  f(t)  (29) 

—oo 

y  en  consecuencia: 

/(f)®8(í)=/(f)  (30) 

Esto  quiere  decir  que  la  convolución  de  una  señal  o 
función  con  un  impulso  unitario  devuelve  como  re¬ 
sultado  la  misma  señal.  Se  trataría  así  el  impulso 
como  el  elemento  neutro  de  la  operación  de  convolu¬ 
ción. 

Pasemos  ahora  a  las  coordenadas  espaciales  x  e  y,  en 
lugar  de  la  variable  lineal  t.  Como  el  espacio  objeto 
es,  en  general,  diferente  al  espacio  imagen,  establez¬ 
camos  para  él  un  sistema  de  coordenadas  r\.  Pues 
bien,  si  consideramos  un  sistema  de  formación  de 
imágenes  X  que  cumpla  las  propiedades  de  lineali- 
dad,  su  respuesta  h{x,y )  ante  la  excitación  producida 
por  un  impulso  unitario  será: 

3ó[5(£,,r|)]  =  h(x,y)  (31) 

Y  si  el  sistema  es  invariante  o  de  parámetros  constan¬ 
tes,  cuando  la  señal  de  entrada  sufra  un  desplaza¬ 
miento  espacial,  la  imagen  de  salida  registrará  el 
mismo  desplazamiento  sin  modificación  de  la  fun¬ 
ción: 

3C  n  -  rj')]  =  h(x  -x',y-  y')  (32) 

Por  la  propiedad  del  impulso  se  puede  expresar  cual¬ 
quier  función  de  entrada  /e(x, y)  como  la  convolución 
de  ella  misma  con  un  impulso  unitario,  en  este  caso 
extendida  la  integral  de  convolución  a  las  dos  dimen¬ 
siones  planimétricas: 
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(33) 


/(^'n)  =  /(^'n)®s(^ri)  = 

00  00 

=  ||  /(w,v)S(£,-w,r|-v) dudv 

"00  "00 


La  función  salida  o  respuesta  del  sistema  formador 
de  imágenes  a  la  señal  de  entrada  es: 

g(x,y)  =  X[/(^,r|)]  = 


=  X 


uu  uu 

v)  8  (H,  -  u,  r|  —  v)  dudv 


(34) 


Y  si  el  sistema  es  lineal  e  invariante: 


X 


v)  8  (H,  -  u,  r|  -  v)  dudv 


00  00 

=  J  J /(w,v)3€[8(^  —  u,r\— vj^dudv  = 

"00  "00 


00  00 

=  ||  /(w,v)-/z(x-w,  y -v^ dudv 

"00  "00 


(35) 


Esta  última  expresión  es  la  de  la  convolución  entre 
la  señal  de  entrada  y  la  función  h(x,y),  que  es  la  res¬ 
puesta  del  sistema  ante  un  impulso  unitario,  de  modo 
que  finalmente  se  puede  escribir: 

g(x,y)=fi(x,y)0/(^,r |)  (36) 


Es  decir,  que  la  respuesta  de  un  sistema  lineal  con¬ 
siste  en  la  convolución  de  la  entrada  con  la  respuesta 
del  sistema  al  impulso  unitario  y  por  lo  tanto,  está 
determinada  por  esta  última.  Puede  conocerse  con 
exactitud  la  respuesta  del  sistema  ante  una  entrada 


de  gran  complejidad  sólo  conociendo  cómo  trans¬ 
forma  el  sistema  una  entrada  unitaria.  En  la  forma¬ 
ción  de  imágenes,  la  función  h(t)  es  la  PSF  vista  an¬ 
teriormente  y  la  imagen  formada  es  en  realidad  la 
convolución  entre  la  imagen  real,  la  correspondiente 
al  espacio  objeto,  y  la  PSF.  Cuanto  mayor  sea  la  dis¬ 
persión  atmosférica,  unida  a  las  distorsiones  causa¬ 
das  por  la  óptica  del  instrumento,  mayor  será  el 
efecto  emborronador  ( blurring )  de  la  PSF  y  peor  será 
la  nitidez  de  la  imagen. 

Se  denomina  función  del  sistema  a  la  transformada 
de  Fourier  de  la  respuesta  del  sistema  al  impulso  uni¬ 
tario: 


//(p,a)  =  3:[/z(x,y)]=  J  jh(x,y)e~,{sa+ay)dxdy  (37) 


La  función  del  sistema  lo  caracteriza  completamente 
y  su  conocimiento  permite  predecir  la  degradación  a 
que  se  verá  sometida  la  imagen.  La  transformada  de 
Fourier  de  una  imagen  permite  definirla  en  el  domi¬ 
nio  de  la  frecuencia  en  lugar  de  hacerlo  en  el  dominio 


Grabación 


Los  satélites  adquieren  líneas  completas  de  imagen 
sucesivamente  y  sin  solución  de  continuidad  mien¬ 
tras  se  desplazan  en  su  órbita  alrededor  de  la  Tierra. 
La  escena  formada  por  un  cierto  paquete  de  líneas  se 
forma  en  el  segmento  terrestre  de  la  misión  para  co¬ 
mercializar  la  información.  Las  celdas  que  constitu¬ 
yen  cada  línea  de  imagen  pueden  ser  adquiridas  si¬ 
multáneamente  o  mediante  barrido  secuencial,  según 
sea  el  diseño  del  sensor.  Incluso  las  líneas  pueden  ser 
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espacial  convencional.  Una  de  las  ventajas  de  operar 
en  el  dominio  de  la  frecuencia  es  que,  en  virtud  del 
teorema  de  la  convolución,  la  operación  de  convolu¬ 
ción  de  dos  funciones  en  el  dominio  espacial  equivale 
al  producto  de  sus  respectivas  transformadas  de  Fou¬ 
rier,  y  multiplicar  imágenes,  aunque  sean  éstas  en  el 
dominio  de  la  frecuencia  es  una  operación  más  sim¬ 
ple  y  permite  mayor  flexibilidad.  Todo  ello  se  verá  en 
mayor  detalle  en  el  apartado  dedicado  a  la  imagen  en 
el  dominio  de  la  frecuencia,  dentro  de  la  lección  de 
tratamiento  digital  de  imágenes. 

Adicionalmente,  si  se  conoce  la  PSF  o  se  tienen  me¬ 
dios  para  estimarla,  operando  en  el  dominio  de  la  fre¬ 
cuencia,  aunque  también  en  el  espacial,  es  posible 
aplicar  funciones  capaces  de  introducir  el  efecto 
opuesto  al  inducido  por  aquella  en  la  imagen.  De  esta 
manera  se  pueden  restaurar  imágenes  degradadas 
por  la  PSF  para  obtener  otras  más  parecidas  a  la 
realidad  aplicando  funciones  restauradoras.  Estos 
procedimientos  también  se  conocen  genéricamente 
como  deconvolución,  haciendo  referencia  al  efecto 
contrario  al  que  causa  la  PSF  sobre  la  imagen,  y  se 
verán  en  el  capítulo  10. 


adquiridas  una  a  una  o  bien  en  bloques  de  varias  de 
ellas. 

El  instrumento  está  dotado  de  una  óptica  y  unos  sis¬ 
temas  electrónicos  que  permiten  dirigir  y  registrar  la 
señal  recibida  en  su  campo  de  visión  en  distintos  pla¬ 
nos  focales  dedicados  a  cada  longitud  de  onda.  En 
cada  uno  de  ellos,  el  dispositivo  se  encarga  de  tradu¬ 
cir  la  irradiancia  espectral  en  niveles  digitales  de  cada 
celda  para  configurar  la  imagen  digital.  Finalmente, 
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esa  imagen  se  graba  a  bordo  del  satélite  y  se  trans¬ 
mite  a  la  estación  de  recepción  terrestre,  bien  direc¬ 
tamente,  bien  a  través  de  constelaciones  de  satélites 
repetidores. 


Ruido 

Se  denomina  ruido  electrónico  a  la  energía  electro¬ 
magnética  que  aparece  en  cualquier  circuito  electró¬ 
nico,  además  de  la  señal  útil  a  procesar.  La  calidad 
de  la  imagen  adquirida  con  cualquier  sistema  forma- 
dor  de  imágenes  está  en  general  afectada  por  un  de¬ 
terminado  nivel  de  ruido,  cuya  existencia  viene  a  per¬ 
turbar  el  procesamiento  e  interpretación  de  la  señal 
original.  Por  extensión,  se  considera  ruido  a  toda  in¬ 
formación  espuria  añadida  por  cualquier  vía  a  la  in¬ 
formación  verdaderamente  útil  que  procede  de  la  su¬ 
perficie  observada  y  que,  en  su  caso,  generaría  una 
reproducción  fidedigna  en  el  espacio  imagen  de  la 
realidad  del  espacio  objeto.  Lo  cierto  es  que  durante 
el  proceso  de  registro  de  la  imagen  se  introduce  inevi¬ 
tablemente  ruido  como  un  factor  aditivo,  y  por  esta 
razón  el  sistema  en  su  conjunto  no  puede  conside¬ 
rarse  exactamente  lineal. 

g(x,y)  =  Grab[f(^  p)]  =  (3g) 

=  h(x,y)®f($,T] |)  +  n(5,ri) 

Mediante  operadores  de  deconvolución  se  puede  res¬ 
taurar  la  imagen,  si  bien  no  completamente,  ya  que 
el  término  de  ruido  es  aditivo  y  no  está  dentro  de  la 
función  de  convolución  degradante  de  la  imagen.  Al 
aplicar  operadores  de  restauración  es  inevitable  per¬ 
der  parte  de  la  señal  útil  mientras  se  busca  eliminar 
solo  el  ruido.  Por  lo  tanto,  en  la  restauración  de  una 
imagen  se  buscará  mantener  un  compromiso  entre 
señal  y  ruido  acorde  con  las  exigencias  del  trabajo. 


Formatos 


La  imagen  multiespectral,  como  se  ha  visto,  consiste 
en  un  pila  de  imágenes  monocanal.  En  los  sistemas 
de  almacenamiento,  discos,  cintas,  memorias  USB, 
etc.  la  imagen  es  un  archivo  de  carácter  lineal  donde, 
si  una  celda  se  almacena  en  un  byte,  el  registro  de  la 
imagen  completa  consiste  en  un  número  determi¬ 
nado  de  bytes,  idealmente  contiguos,  igual  al  número 
de  celdas  contenido  en  el  conjunto  de  todas  las  ban¬ 
das.  Existen  tres  formatos  de  grabación  para  este  tipo 
de  imágenes:  bandas  secuenciales  (BSQ),  bandas  in¬ 
tercaladas  por  línea  (BIL)  y  bandas  intercaladas  por 
píxel  (BIP).  En  el  primero  de  ellos  se  escribe  la  ima¬ 
gen  registrando  todas  las  celdas  de  una  en  una  por 
líneas  consecutivas  y  una  banda  tras  otra.  Cuando  el 
almacenamiento  se  realizaba  en  los  inicios  del  trata¬ 
miento  digital  sobre  cintas  magnéticas,  estudiar  un 
determinado  píxel  multiespectral  requería  la  bús¬ 
queda  del  byte  en  cuestión  a  lo  largo  de  toda  la  cinta, 
ya  que  una  banda  no  era  grabada  hasta  no  haber  ter¬ 
minado  completamente  la  anterior.  Lógicamente  el 
tiempo  de  proceso  se  dilataba  considerablemente. 
Para  agilizar  estos  procesos  se  recurrió  al  tercero  de 
los  sistemas  (BIP)  y  comenzaron  a  almacenarse  los 
niveles  digitales  de  cada  celda  en  todas  las  bandas  de 
la  imagen  en  las  posiciones  sucesivas  del  archivo  in¬ 
formático.  Concluidos  los  registros  de  los  n  valores 
de  la  celda  en  las  n  bandas,  se  pasaba  a  la  siguiente 
celda,  y  así  sucesivamente.  Si  con  este  procedimiento, 
el  acceso  a  un  píxel  concreto  era  más  rápido,  en  cam¬ 
bio  la  lectura  de  una  banda  retrasaba  notablemente 
el  proceso.  Una  solución  de  compromiso  consistió  en 
grabar  una  línea  completa  de  una  banda,  pasar  a  con¬ 
tinuación  a  grabar  la  misma  línea  de  la  segunda 
banda,  y  así  sucesivamente  hasta  completar  todas  las 
bandas.  Terminado  el  registro  de  la  primera  línea 
multiespectral  se  pasaba  a  la  segunda  línea  y  se  con¬ 
tinuaba  el  proceso  hasta  terminar  todas  las  líneas  de 
la  imagen  (BIL).  En  la  actualidad  coexisten  los  tres 


Formatos  de  grabación  de  la  imagen  digital. 


Carlos  Pinilla 


sistemas  y  los  paquetes  de  tratamiento  digital  de  imá¬ 
genes  permiten  cómodamente  la  codificación  de  cual¬ 
quiera  de  ellos  y  la  traducción  de  uno  a  otro. 


Calibración 


En  determinados  estudios  basados  en  la  teledetec¬ 
ción  no  basta  con  medir  las  diferencias  en  la  res¬ 
puesta  espectral  de  las  diferentes  zonas  de  la  imagen, 
sino  que  es  necesario  conocer  de  modo  absoluto  el 
comportamiento  físico  de  la  cubierta. 

Para  ello  es  necesario  encontrar  la  función  genérica 
de  transformación  que  existe  entre  la  variable  física 
objeto  de  estudio  y  el  ND  asignado  finalmente  por  el 
sistema  de  adquisición.  Una  vez  conocida  esa  fun¬ 
ción,  la  aplicación  de  su  inversa  a  los  ND  de  la  imagen 
permite  recuperar  la  variable  física  y  poder  así  carac¬ 
terizar  la  cubierta  estudiada 

Calibrar  una  imagen  consiste  es  transformar  sus  ni¬ 
veles  digitales  en  variables  que  tengan  un  significado 
físico  como  son  la  radiancia,  la  reflectancia,  la  tempe¬ 
ratura  de  la  superficie,  etc. 

La  relación  entre  ND  y  la  magnitud  física  habitual¬ 
mente  implica  dos  etapas: 

1.  Convertir  los  ND  en  valores  de  radiancia  espectral 
y,  tras  ello,  en  reflectancia  espectral  de  la  cubierta: 
La  relación  entre  unos  y  otros  está  gobernada  por 
funciones  matemáticas  simples,  normalmente  li¬ 
neales,  las  llamadas  funciones  de  calibración. 

2.  Establecer  una  relación  entre  la  reflectancia  ha¬ 
llada  y  la  magnitud  física  buscada.  Este  paso  se 
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basa  en  dependencias  estadísticas,  según  la  corre¬ 
lación  que  exista  entre  aquélla  y  ésta 

Función  de  conversión 

Los  sensores  integrados  en  el  instrumento  de  adqui¬ 
sición  convierten  la  irradiancia  recibida  en  una  señal 
eléctrica  variable.  El  voltaje  producido  es  proporcio¬ 
nal  a  la  energía  incidente. 

Dado  que  se  conocen  los  parámetros  orbitales  del  sa¬ 
télite,  y  con  ello  su  situación  con  respecto  a  la  super¬ 
ficie  observada,  particularmente  su  distancia  y  el  án¬ 
gulo  de  observación,  es  inmediato  deducir  la  radian¬ 
cia  que  procede  del  terreno  a  partir  de  la  irradiancia 
sobre  el  plano  focal  del  instrumento.  En  consecuen¬ 
cia,  se  puede  decir  que  el  sensor  convierte  la  radian¬ 
cia  en  una  señal  analógica,  materializada  por  una  co¬ 
rriente  eléctrica. 

Otro  elemento  constitutivo  del  instrumento  es  un 
convertidor  analógico-digital  que  transforma  esa  co¬ 
rriente  en  niveles  digitales  discretos,  en  función  del 
voltaje  de  la  primera. 

Existen  diversos  tipos  de  sensores,  cuya  descripción 
queda  fuera  del  alcance  de  este  manual.  Los  más  im¬ 
portantes  son  de  tipo  CCD  (Charge-Coupled  Device, 
dispositivo  de  carga  acoplada)  y  los  de  tipo  CMOS 
( Complementar y  Metal-Oxide  Semiconductor,  semi¬ 
conductor  complementario  de  óxido  metálico). 
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Basándose  ambos  en  el  efecto  fotoeléctrico,  el  detec¬ 
tor  contiene  tantos  fotodiodos  o  fotosites  como  píxe- 
les,  los  cuales  producen  una  corriente  eléctrica  pro¬ 
porcional  al  número  de  fotones  incidentes  y,  en  defi¬ 
nitiva,  a  la  intensidad  de  luz  recibida. 

En  los  sensores  CCD,  esa  corriente  se  transfiere  a  un 
circuito  integrado  que  contiene  una  serie  de  conden¬ 
sadores  acoplados,  cada  uno  de  los  cuales  se  carga 
en  función  de  la  corriente  recibida.  En  electrónica  in¬ 
dependiente,  la  señal  se  amplifica,  y  en  el  convertidor 
analógico-digital  la  carga  inicial  es  convertida  en  ni¬ 
veles  digitales  que  son  finalmente  almacenados  en 
forma  de  imagen  digital. 

Los  sensores  CMOS  son  en  realidad  sensores  de  pí- 
xeles  activos  basados  en  la  tecnología  CMOS.  A  dife¬ 
rencia  del  CCD  incorpora  un  amplificador  de  la  señal 
eléctrica  en  cada  fotodiodo  y  suele  incluir  el  conver- 
sor  analógico-digital  en  el  propio  chip,  prescindiendo 
así  de  parte  de  la  electrónica  externa. 

Los  sensores,  con  independencia  de  la  familia  a  la 
que  pertenezcan,  transforman  en  una  señal  de  salida 
Q(V),  materializada  por  un  voltaje  V,  la  radiancia  es¬ 
pectral  de  la  superficie  L\,  medida  en  W-tn2-sr\  Esta 
correspondencia  es  prácticamente  lineal  en  el  inter¬ 
valo  espectral  de  cada  una  de  las  bandas  exploradas 
por  el  sensor. 

Q  =  k\dXLxhdX  +  0  (1) 

Donde  k  es  un  factor  de  proporcionalidad,  n  es  la 
respuesta  espectral  del  instrumento  y  O  es  una  cons¬ 
tante  aditiva.  De  forma  simplificada  puede  escribirse: 

Q  =  G1Lx+Ol  (2) 


de  donde  se  deduce  que: 


Q-0 


G 


Es  decir: 


e  =  /,(4);  4  =  /r'(e)  e» 

El  voltaje  es  función  de  la  radiancia  y  viceversa,  y, 
como  puede  verse,  la  relación  entre  la  señal  de  salida 
y  la  radiancia  tiene  la  forma  de  una  recta. 

De  otra  parte,  la  relación  entre  la  radiancia  y  el  nivel 
digital  es  también  lineal: 


ND  —  G2LX  +  O 2 


(4) 


y 

r  _nd-o2 
h=  g2 

por  lo  que: 

ND  =  f2(Lx );  Lx  =  f2~l(ND )  (5) 

Es  decir,  también  puede  ponerse  en  relación  directa¬ 
mente  mediante  funciones  de  calibración  la  radiancia 
del  terreno  y  el  nivel  digital  de  cada  una  de  las  celdas 
de  una  imagen.  En  la  práctica,  los  sensores  se  cali¬ 
bran  exponiéndolos,  al  menos,  a  dos  niveles  conoci¬ 
dos  de  radiancia  y  calculando  los  coeficientes  de  ga¬ 
nancia  G  ( gain )  y  sesgo  O  [offset).  Lo  habitual  es  de¬ 
ducir  la  función  de  calibración  por  ajuste  a  la  recta 
de  regresión  del  conjunto  de  numerosas  lecturas. 
Estas  operaciones  se  realizan  en  tierra  durante  la 
fase  de  prelanzamiento  del  satélite,  antes  de  su 


O _ 

L 

Linealidad  del  voltaje  con  respecto  a  la  radiancia. 


Radiancia 


Funciones  de  conversión  de  radiancia  en  nivel  digital. 

Q 


Q  =  GL+0 
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Ajuste  por  regresión  de  los  parámetros  de  calibración. 


Calibración 
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puesta  en  órbita.  Pero  una  vez  orbitando,  vuelve  a 
recalibrarse  periódicamente  para  ajustar  los 
parámetros  de  ganancia  y  sesgo,  ya  que  por  el 
envejecimiento  de  los  materiales  suelen  variar, 
aunque  mínimamente.  Esta  recalibración  se  lleva  a 
cabo  realizando  medidas  simultáneas  de  radiancia 
desde  el  satélite  y  desde  tierra,  mediante 
espectrorradiómetros,  sobre  blancos  conocidos. 


Distintos  modelos  de  espectrorradiómetros. 


Radiometría  de  campo 

La  radiometría  consiste  en  medir  la  radiación  emitida 
o  reflejada  desde  una  determinada  superficie  objeto 
de  estudio,  mediante  espectrorradiómetros,  sometién¬ 
dola  a  distintas  condiciones  de  observación.  Las 
campañas  de  radiometría  son  útiles  y  necesarias  para 
caracterizar  desde  el  punto  de  vista  raiométrico  la 
verdad-terreno  ( ground-truth )  con  distintos  objetivos 
como  son  la  calibración  de  los  intrumentos  de 
adquisición,  la  validación  de  los  datos  adquiridos  por 
los  satélites  y  la  interpretación  de  los  resultados. 

El  radiómetro  es  un  instrumento  destinado  a  realizar 
medidas  en  campo  o  en  laboratorio,  de  lo  cual  depen¬ 
derá  su  tamaño  y  complejidad,  capaz  de  medir  la  ra¬ 
diancia  espectral  en  un  cierto  intervalo  de  longitudes 
de  onda.  Suele  ser  un  instrumento  hiperespectral, 
esto  es,  realiza  medidas  de  radiancia  en  intervalos 
muy  estrechos  y  contiguos  de  longitudes  de  onda,  fre¬ 
cuentemente  de  entre  1  y  5  nm. 

Mediante  la  radiometría  de  campo  es  habitual  reali¬ 
zar  medidas  bidireccionales,  con  el  fin  de  caracterizar 
el  hábito  reflexivo  de  la  superficie,  su  eventual  carác¬ 
ter  lambertiano,  deducir  la  función  de  distribución  de 
reflectancia  bidireccional,  etc.  Para  ello  se  mueve  el 
radiómetro  haciéndolo  girar  dentro  de  un  plano  ver¬ 
tical,  que  a  su  vez  también  puede  girar  alrededor  del 
eje  normal  a  la  superficie.  De  esta  manera,  mante¬ 
niéndose  la  fuente  de  iluminación  fija,  el  sol,  se  mide 


la  radiancia  en  todas  las  direcciones  del  hemisferio 
superior  a  la  superficie. 

También  es  frecuente  realizar  medidas  manteniendo 
fijo  el  espectrorradiómetro  y  dejando  que  el  sol 
realice  su  órbita  diurna,  a  fin  de  comprobar  la  res¬ 
puesta  de  la  superficie  bajo  distintas  condiciones  de 
azimut  y  elevación  solar. 

Calibración  de  sensores  orbitales 

Genéricamente,  la  función  de  calibración  de  un  sen¬ 
sor  formador  de  imágenes  situado  en  órbita  tiene  la 
forma: 

L,„{k)  =  a,  (k)  +  a,  (k)ND(k)  (6) 

en  cuya  expresión  Lsen  es  la  radiancia  medida  por  el 
sensor,  ao(k)  y  ai(k)  son  los  coeficientes  de  calibra¬ 
ción  del  sensor  y  ND(k)  es  el  nivel  digital,  todo  ello 
referido  a  una  banda  genérica  k  de  la  imagen. 

En  el  espectro  óptico,  la  emisión  propia  de  la  super¬ 
ficie  terrestre  es  despreciable  y,  por  ello,  se  considera 
que  la  radiación  procedente  de  la  tierra  en  ese  inter¬ 
valo  espectral  se  debe  a  la  reflexión.  Los  valores  de 
radiación  registrados  dependen  entonces  de  cuatro 
factores: 

-i  La  irradiancia  solar,  que  depende  a  su  vez  de 
tanto  de  la  altura  solar  como  de  la  distancia 
Tierra-Sol  en  el  momento  de  la  adquisición  de 
la  imagen. 

-i  La  reflectancia  de  la  superficie,  que  es  un  pa¬ 
rámetro  invariable,  siempre  que  la  naturaleza 
o  la  cobertura  de  la  superficie  se  mantenga 
en  el  tiempo. 

-i  La  transmitancia  atmosférica  en  las  longitu¬ 
des  de  onda  consideradas,  que  depende  de  la 
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composición  química  de  los  gases  atmosféri¬ 
cos  y  del  contenido  en  aerosoles. 

-i  La  radiancia  difusa  atmosférica,  que  depende 
de  los  aerosoles  en  suspensión. 

Puede  decirse  que,  en  su  tránsito  a  través  de  la  at¬ 
mósfera,  parte  de  la  energía  solar  incidente  sobre  el 
borde  exterior  de  la  atmósfera  se  reduce,  debido  a 
fenómenos  de  absorción,  hasta  llegar  a  la  superficie, 
donde  es  reflejada  en  virtud  de  la  reflectancia  inhe¬ 
rente  al  terreno  y  redirigida  en  numerosas  direccio¬ 
nes,  entre  ellas  también  la  del  sensor.  De  nuevo,  al 
atravesar  la  atmósfera,  parte  de  la  energía  es  absor¬ 
bida,  pero  el  fenómeno  se  ve  compensado  parcial¬ 
mente  por  la  adición  de  la  radiancia  difusa  atmosfé¬ 
rica  que  ensucia  la  señal  y  termina  penetrando  en  el 
campo  de  visión  del  instrumento. 

La  radiancia  medida  por  el  sensor  es: 


 E0(k) eos Q 
K% 


En  donde  Eo(k )  es  la  irradiancia  espectral  solar  extra- 
terrestre  en  la  banda  k,  0  es  el  ángulo  cenital  solar  y 
K  es  un  factor  corrector  de  la  distancia  Tierra-Sol  que 
puede  obtenerse  de  un  anuario  astronómico.  Combi¬ 
nando  las  dos  expresiones  anteriores: 


p(jfc)  =  Kn  Lsen^  CJk) 
E0(k)x(k)cosQ 


(10) 


Si  se  corrige  el  efecto  atmosférico,  Lam(k)= 0  y  t(ft)=l, 
por  lo  que: 


p(&)  =  Kn 


4.W 

E0(k)  cos0 


(11) 


m  =  Lsen(k)-Latm(k)  =  p(k)Lsue(k)x(k) 


(7) 


donde  Lsen(k)  es  la  radiancia  espectral  medida  en  el 
sensor,  Lsue{k)  es  la  radiancia  espectral  solar  recibida 
sobre  el  suelo,  p(k)  es  la  reflectancia  espectral  de  la 
cubierta,  x{k)  es  la  transmitancia  espectral  de  la  at¬ 
mósfera  y  Latm{k)  es  la  radiancia  espectral  aportada 
por  la  atmósfera. 


Y  considerando  (6): 

a0  (k)  +  ax  ( k)ND(k ) 


p  (k)  =  Kn- 


E0(k)  cos0 


(12) 


expresión  que  puede  ponerse  de  la  forma: 
p  (k)  =  b0(k)  +  bl(k)ND(k) 


(13) 


Despejando: 

Lm(k)-Lmm(k)  <8» 

PW-  L,Jk)z(k) 

Si  se  considera  la  superficie  como  si  fuera  lamber- 
tiana,  para  un  determinado  día  del  año,  la  radiancia 
incidente  sobre  el  suelo  es: 


y  que  indica  que,  de  forma  análoga  a  la  radiancia, 
también  la  reflectancia  muestra  una  relación  lineal 
con  el  nivel  digital  de  cada  celda  en  la  imagen,  a  tra¬ 
vés,  en  este  caso,  de  los  coeficientes  de  calibración 
bo(k )  y  b\(k )  de  cada  banda.  Al  igual  que  los  coefi¬ 
cientes  de  calibración  de  radiancia,  los  de  reflectancia 
son  proporcionados  por  el  constructor  del  sensor  y 
se  actualizan  periódicamente. 


Transferencia  radiativa  a  través  de  la  atmósfera. 
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